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Flat spectrum radio quasars (FSRQs) are bright active galactic nuclei surrounded by gas clouds within a UV- 
visible intense radiation field that form the so-called broad line region (BLR). These objects emit relativistic 
jets from a region close to the central supermassive black hole and through the BLR. The Fermi-haige Area 
Telescope (Fermi-LAT) is sensitive to gamma-ray photons from ~30 MeV to more than 300 GeV. We have 
performed spectral analysis of bright FSRQs in a 5.5 year (2008-2014) data sample collected by Fermi-LAT, 
using the new Pass 8 event selection and instrument response function. Also, our study of flaring episodes in a 
limited time range brings interesting results while compared to the full 5.5 year data samples. 


1. Modelling the BLR radiation fieid 


FSRQs constitute a class of active galactic nuclei 
(AGNs) with a dense BLR in which gamma rays 
with energy > 10 GeV are absorbed due to electron- 
positron pair creation, if produced deep inside the 
BLR. Indeed, BLR is expected to be denser in FS¬ 
RQs compared to the BL Lac class. 

Operating since 2008, the Fermi satellite has 
amassed more than 6 years of data, continuously sur¬ 
veying the whole sky Q. The sensitivity of Fermi- 
LAT is ideal for the study of the gamma-ray absorp¬ 
tion inside FSRQs in the 100 MeV-300 GeV range. 
From constraints on gamma-ray absorption we may 
infer limits on the location of the gamma-ray emis¬ 
sion region in the FSRQ jets. 

We expect > 10 GeV photons of FSRQs to un¬ 
dergo absorption in the BLR, where the target photon 
with energy e is a UV photon from the BLR radiation 
field. As most of these photons are expected to come 
from the emission lines, we use a model that includes 
the 6 strongest lines (NV, Lya, OVI Ly/3, GUI NIII, 
NeVIII OIV, Hell Lya) between ~10 to 41 eV [H. 

We model these lines using a Breit-Wigner distri¬ 
bution, given by: 


BW{€) 
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where Ui and lUi are the number density and width, 
respectively, for a given line i. 

Under the commonly used relations Lblr = 
0-1 ^disc and Rblr = VLdisc i 0, where Lblr 
is the luminosity of the BLR, Rblr its radius, and 
Ldisc the luminosity of the accretion disc. The pho¬ 
ton density of the radiation field for each line i can 
be written: 
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The opacity is derived from and is expressed 
as a function of E and z: 
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where tq is the classical electron radius, m the elec¬ 
tron mass. 

The Ty-y opacity in the BLR is then calculated as 
following: 


T^-f[E, z) = ax Rblr x W 

assuming the gamma rays are produced within 
a X Rblr, where Rblr is the outer radius of BLR, 
and a < 1. Since this absorption happens at some dis¬ 
tance from the supermassive black hole, this corrective 
factor that we called “a” represents the fraction of the 
BLR responsible for the absorption. 

In Table U are displayed the line properties of the 
average spectrum of quasars we used in our model, 
as they were given in [I^. Since the He H Lya line 
has quite large uncertainties, we arbitrary fixed its 
EW and relative flux to be equal to the ones of N V 
(uncertainty represented by (*)). 

Very high energy gamma rays travelling from 
far distances undergo absorption in the extragalac- 
tic background light (“EBL”, mainly composed of 
infrared-UV radiation). This absorption is to be con¬ 
sidered above 10 GeV and has been implemented in 
our studies, from the model presented in Q. 

Evidence of absorption in the BLR for some FSRQs 
have been reported in 
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Opacity in the BLR using the Breit-Wigner line model 



Figure 1: Opacity versus Energy, for 3C 454.3. 

The opacity sum on the 6 lines is drawn in plain black. 


Table I Properties of 5 main lines of the average 
spectrum of quasars as compiled in 0, and He II Lya 
as we defined it for this study. 


Line 

e (eV) 

EW (eV) 

Relative flux 

NV 

10.0 

0.16 

0.22 

Lya 

10.2 

0.71 

1.00 

OVI Lyp 

12.04 

0.19 

0.191 

cm NIII 

12.65 

0.09 

0.081 

NeVIII OIV 

15.90 

0.08 

0.047 

Hell Lya 

40.81 

0.16* 

0.22* 


2. Data processing and model fitting 

We have analysed data of 7 bright gamma-ray 
FSRQs. Plots of the spectral energy distributions 
(SEDs) under the label “5.5 years” have been pro¬ 
cessed from 4 August 2008 till 30 April 2014. The 
sources we present in this paper are listed in Table |TT1 
LI I and BII being respectively the Galactic latitude 
and longitude in decimal degrees. 

Data were processed using the Pass 8 data repre¬ 
sentation (P8_SOURCE_V4), and the Science Tools 
version v9-34-01. Signal is reconstructed from each 
source using the unbinned likelihood tool^, applied to 
LAT data in the 0.1-300 GeV energy range, within a 
region of interest (ROI) of 10° radius. A source region 
extended to an additional 10° annulus accounted for 
all the point sources of the Second Fermi-LAT source 
catalog [9j , and for the Galactic diffuse emission (tem- 
plate_4years_P8_V2_scaled) and the isotropic diffuse 
emission (isotropic_source_4years_P8V3). 

We computed the SEDs for all the sources of the 
selected sample with the Pass 8 data representation. 
Additionally, for the two brightest objects of our 


^ http: / / fermi.gsfc.nasa.gov/ssc / data/analysis / 


3C454.3 /04 Aug 2008 - 20 Jan 2014 / Pass 7 vs Pass 8 



Energy (MeV) 

Figure 2: Comparison of the SED of 3C 454.3, produced 
with the Pass 7 (magenta) and Pass 8 (green) data sets. 


FSRQ sample, i.e. 3G 454.3 and PKS 1510-089, we 
also computed the SEDs with the PASS 7 reprocessed 
dataset (P7REP_SOURGE_V15) and verified the con¬ 
sistency of the results with respect to the PASS 8 ones. 
Although some bin-to-bin fluctuations appear due to 
energy wise event migrations, the two SEDs (Pass 7 
and Pass 8) for both the sources are compatible (as 
shown in Figure [2] for 3G 454.3). 

A first set of fits was performed from 100 MeV 
till the highest energy data point (excluding up¬ 
per limits), using a log-parabola (LP: dN{E)/dE = 
No {E/Eo)-°‘-f^ iog(E/Eo))^ gxed at 

297.6 MeV, and where “log” is the natural loga¬ 
rithm), a broken power law (BPL: dN{E)/dE = 
No {E/Eb)-^% with i = I a E < Et 
and i = 2 if A > Eb), and a power law 
with an exponential cutoff (PLEC: dN{E)/dE = 
No exp{—E/Ec), with Ep kept fixed 

at 412.7 MeV). Other sets of fits were performed by 
addiM exponential factors to model the EBL absorp¬ 
tion p and the opacity of the BLR. 

As the fits presented in the Sections [3] and 2] are bin¬ 
ning dependant, the values of the fit parameters vary 
from one choice of binning to another. Narrow data 
binning could held spurious fluctuations, while wide 
data-binning could hide features. In order to estimate 
this systematic effect, we do the following: a first LP 
fit is performed on the SED, while keeping all param¬ 
eters fixed to the values obtained by the unbinned 
likelihood analysis, and a is returned. A sec¬ 

ond LP fit is performed with No, a and /3 kept free, 


eConf C141020.1 




















































5th Fermi Symposium : Nagoya, Japan : 20-24 Oct, 2014 


3 


Table II Characteristics of the 7 FSRQs used in this paper. Sources are ordered by decreasing flux in the 1-100 GeV 
energy range. Values of the luminosity in the broad line region {Lblr) are taken or derived from [ia,[ii- 


Name 

RA 

(J2000.0) 

Dec 

(J2000.0) 

LII 

(deg.) 

BII 

(deg.) 

Redsliift 

z 

Flux 2FGL 

1-100 GeV 

Flux IFHL 

10-100 GeV 

Photon index 

IFGL 

Lblr 

{10'^erg.s~^) 

3C 454.3 

22 53 57.7 

+16 08 53.1 

86.11 

-38.19 

0.859 

9.65e-8 

1.35e-9 

2.46619 

33.00 

PKS 1510-08 

15 12 50.5 

-09 06 00.9 

351.28 

40.13 

0.360 

4.06e-8 

7.35e-10 

2.40756 

5.62 

4C +21.35 

12 24 54.5 

+21 22 46.9 

255.08 

81.65 

0.434 

3.54e-8 

7.43e-10 

2.54717 

15.80 

3C 279 

12 56 11.0 

-05 47 20.1 

305.1 

57.06 

0.536 

2.56e-8 

5.37e-10 

2.32061 

3.10 

PKS 0454-234 

04 57 03.1 

-23 24 52.0 

223.7 

-34.9 

1.003 

2.27e-8 

2.99e-10 

2.20649 

3.70 

B2 1520+31 

15 22 09.8 

+31 44 14.3 

50.16 

57.02 

1.484 

1.76e-8 

4.27e-10 

2.42125 

8.00 

PKS (B) 1424-418 

14 27 56.2 

-42 06 18.6 

321.44 

17.26 

1.522 

1.47e-8 

2.9e-10 

2.31004 

8.91 


Eb kept fixed to the value returned by the unbinned 
likelihood analysis. We compare the fits to validate 
that the binned fit is compatible with the unbinned 
fit, though results differ. These systematics could be 
overcome by implementing an unbinned analysis for 
all fitted models in future. 

Our modelling study is then done by fitting the 
SEDs using the LP, BPL and PLEC function that all 
now include EBL. These fits will be reported in the 
Sections [3] and 01 and compared to the fits that in¬ 
clude both EBL and BLR absorption. The latter fits 
are written as LPt, BPLt, and PLECt. 

The observed spectrum Fobs{E) will then be ulti¬ 
mately written: 

Fobs{E) = e-“ F,^t{E), (5) 

where Fint{E) is the LP, BPL or PLEC fitting func¬ 
tion. Parameter a is kept free in the [10“^, 1] range to 
account for the fraction of radius of the BLR in which 
gamma rays may be absorbed. 

The fitting procedure using the absorption models 
is implemented by interpolation of both the tebl and 
Tjy graphs. While comparing each “EBL -|- BLR ab¬ 
sorption” fit (LPr, BPLr or PLECr) with the “EBL 
-I- no BLR absorption” fit (LP, BPL or PLEC), and 
if both fits have a x^/nd/ ^ 1, we obtain a p-value 
which indicates the discrepancy between the fit with 
model and the fit without model, for a given function. 


3. Results on the 5.5 years of data 

We present now the results on the 7 bright FSRQs 
during the 5.5 years period previously defined in Sec¬ 
tional Data under the label “5.5 years” are processed 
from 4 August 2008 till 30 April 2014. In Figure 01 
are presented the SEDs of 3C 454.3, PKS 1510-089, 
TXS 1520-k31 and PKS 1424-41. Some of the fits are 
not visible, mainly the BPL and PLEC. The reason 
is that they are overwritten by the BPLr and PLECr 
functions, for which the parameter a is very small. 


We would consider having evidence for absorption 
in the BLR if we get all of the following: 

• at least one good quality fit among one of the fits 
with EBL-I-BLR absorption (LPr, BPLr and 
PLECr); 

• parameter “a” with a relatively small error bar; 

• a small p-value or a bad fit of the corresponding 
function with only EBL absorption (LP, BPL or 
PLEC). 

In Table Uni are displayed the fit parameters of the 
7 sources, along with the p-values used to compare 
the models (with versus without absorption). In blue 
bold face are the parameters that suggest a possible 
BLR absorption, as some of the above conditions are 
partially met for 3C 454.3 (with a p-value of 6.2 x 10“^ 
/ 3.9 (T C.L.), and for TXS 1520-1-31 (with a p-value of 
9.3 X 10“^ / 2.6 a C.L.). We have no hint of absorption 
for the other sources we studied. 


4. Results on high state/flaring episodes 

Data were analysed during flaring/high state peri¬ 
ods for the following sources: 

• 3C 454.3 (high state and giant flare) during 02 
Nov-05 Dec 2010 (MJD 55502.5-55535.5) [H; 

• PKS 1510-089 during 19 Feb-04 Apr 2012 (MJD 
55976.0-56021.0), along with MAGIC data from 
the same period (MAGIC data points taken 
from 0); 

• PKS 1424-41 during 30 Sep 2012-27 Jul 2013 
(MJD 56200.0-56500.0), as a combined series of 
4 successive radio flares). 

During these outburst episodes, the gamma-ray 
emission region can have a different location compared 
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3C 454.3 (Pass 8 - 5.5 years) 


PKS 1510-089 (Pass 8 - 5.5 years) 




TXS 1520+31 (Pass 8 - 5.5 years) 


PKS 1424-41 (Pass 8 - 5.5 years) 




Figure 3: SEDs of 3C 454.3, PKS 1510-089, TXS 1520+31 and PKS 1424-41. BPL and PLEC fits are often hidden 
beneath BPLt and PLECt. 


to the quiescent state. In this section we present the 
results obtained on these three high states (Figure 0] 
and Table HVl) . Under the same criteria than the ones 
used in Section [3l we still find no evidence of BLR ab¬ 
sorption, though we still have a hint of it for 3C 454.3 
with a p-value of around 2.2 x 10“^ for the discrep¬ 
ancy between the BPL and BPLr fits (significance of 
about 3 a). Due to the unusual shape of the SED of 
PKS 1424-41 during this series of 4 flares, all fits have 
a large value. 


Though we dispose of less photons in the data anal¬ 
ysis of flaring episodes, during strong and long flares, 
it could be possible to constrain the location of the 
gamma-ray production region if it is deep enough in¬ 
side the BLR. 
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3C 454.3 (Pass 8 - Nov-Dec2010) 


PKS 1510-089 (Pass 8 - Feb-Apr 2012) 


PKS 1424-41 (Pass 8 - Sep-Jul 2013) 





Figure 4: SEDs of 3C 454.3, PKS 1510-089 and PKS 1424-41 during the outburst periods. PLEC fits are hidden 
beneath PLECt for 3C 454.3 and PKS 1424-41. The flare of PKS 1510-089 was studied along with MAGIC data above 
~90 GeV and the combined LAT-MAGIC SED was fitted only with LP and LPr function (LP hidden beneath LPr). 


5. Conclusions and perspectives 

We find that the gamma-gamma absorption in the 
BLR is not significant enough to claim discovery for 
the models of BLR and spectral functions we have in¬ 
vestigated. There are hints of absorption in case of 3C 
454.3 and TXS 1520-1-31 with significance of the order 
of 3(7. An implication of our results could be that the 
gamma-ray emission zone in FSRQs might be located 
outside or at the outer edge of the BLR. However, fur¬ 
ther investigation on binning effects on the SED fits 
are required. Future work is also expected to improve 
the modelling of the BLR. 
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Parameter 

Model 


3C 454.3 

PKS 1510-089 

4C 4-21.35 

3C 279 

PKS 0454-234 

TXS 1520+31 

PKS1424-24 


LP.unbmned 

OL 

2.248 ± 0.005 

2.294 ± 0.008 

2.240 ± 0.009 

2.254 ± 0.010 

2.067 ± 0.013 

2.273 ± 0.013 

2.090 ± 0.008 


LP.unbmned 

P 

0.0897 ± 0.0024 

0.0513 ± 0.0036 

0.0445 ± 0.0044 

0.0588 ± 0.0053 

0.0793 ± 0.0053 

0.0689 ± 0.0059 

0.0753 ± 0.0039 


PL 

Eb 

298 ± 4.4 

258 ± 6.4 

317 ± 10.8 

284 ± 10.2 

358 ± 16.9 

281 ± 11.2 

620 ± 23.8 

function 

LP 

Oi 

2.240 ± 0.008 

2.292 ± 0.030 

2.274 ± 0.025 

2.266 ± 0.027 

2.070 ± 0.019 

2.274 ± 0.026 

2.077 ± 0.015 

(without BLR 

LP 

p 

0.0980 ± 0.0039 

0.0511 ± 0.0093 

0.0404 ± 0.0092 

0.0539 ± 0.0099 

0.0761 ± 0.0087 

0.0687 ± 0.0110 

0.0713 ± 0.0088 

absorption 

BPL 

Pi 

2.242 ± 0.345 

2.363 ± 0.030 

2.280 ± 0.028 

2.332 ± 0.029 

2.138 ± 0.016 

2.317 ± 0.033 

2.066 ± 0.010 

model) 

BPL 

r2 

2.856 ± 2.108 

2.673 ± 0.075 

2.521 ± 0.057 

2.690 ± 0.098 

2.681 ± 0.064 

2.791 ± 0.121 

2.564 ± 0.087 


BPL 

Eb 

1307 ± 5211.8 

2387 ± 971.5 

1439 ± 573.3 

2738 ± 979.1 

3062 ± 0.3 

2036 ± 674.1 

2947 ± 573.6 


PLEC 

rpLEC 

2.241 ± 0.016 

2.386 ± 0.034 

2.340 ± 0.035 

2.324 ± 0.032 

2.088 ± 0.024 

2.308 ± 0.038 

2.028 ± 0.013 


PLEC 

Ea 

8250 ± 721.2 

40228 ± 13971.0 

52905 ± 28547.3 

26932 ± 8335.2 

14581 ± 2929.5 

14378 ± 3950.2 

16459 ± 2097.1 


LPr 

OL 

2.243 ± 0.008 

2.294 ± 0.033 

2.274 ± 0.025 

2.268 ± 0.029 

2.082 ± 0.017 

2.289 ± 0.033 

2.079 ± 0.013 


LPr 

P 

0.0876 ± 0.0051 

0.0498 ± 0.0117 

0.0404 ± 0.0092 

0.0526 ± 0.0130 

0.0590 ± 0.0103 

0.0506 ± 0.0162 

0.0542 ± 0.0142 

function 

BPLr 

Pi 

2.250 ± 0.012 

2.312 ± 0.037 

2.251 ± 0.031 

2.301 ± 0.031 

2.006 ± 0.025 

2.272 ± 0.045 

2.052 ± 0.011 

(with BLR 

BPLr 

P2 

2.738 ± 0.049 

2.583 ± 0.056 

2.507 ± 0.048 

2.590 ± 0.109 

2.283 ± 0.034 

2.518 ± 0.069 

2.403 ± 0.066 

absorption 

BPLr 

Eb 

1225 ± 142.9 

1152 ± 372.0 

1062 ± 379.5 

1616 ± 743.9 

544 ± 101.1 

769 ± 315.1 

2052 ± 1.2 

model) 

PLECt 

LpLEC 

2.241 ± 0.016 

2.386 ± 0.034 

2.341 ± 0.035 

2.324 ± 0.032 

2.088 ± 0.024 

2.308 ± 0.038 

2.028 ± 0.013 


PLECt 

Ec 

8251 ± 721.4 

40239 ± 13864.7 

52946 ± 27502.4 

26936 ± 8309.2 

14582 ± 2917.7 

14381 ± 3982.3 

16463 ± 2088.6 


LPr 


0.00516 ± 0.00185 

0.00147 ± 0.01118 

0.00001 ± 0.00139 

0.00232 ± 0.02143 

0.02140 ± 0.01007 

0.01500 ± 0.00943 

0.00745 ± 0.00515 

a 

BPLr 


0.00666 ± 0.00195 

0.00785 ± 0.00735 

0.00001 ± 0.00167 

0.00919 ± 0.01496 

0.03641 ± 0.00891 

0.02002 ± 0.00881 

0.00558 ± 0.00488 


PLECt 


0.00001 ± 0.00055 

0.00001 ± 0.00529 

0.00001 ± 0.00151 

0.00001 ± 0.00682 

0.00001 ± 0.00606 

0.00001 ± 0.01584 

0.00001 ± 0.00248 


LP_iin5mned 


36.419 (30) 

7.454 (17) 

18.821 (16) 

5.288 (16) 

14.783 (16) 

8.988 (15) 

36.173 (16) 

X^{n-df) 

LP 


29.316 (27) 

7.278 (14) 

13.009 (13) 

4.609 (13) 

11.516 (13) 

6.448 (12) 

22.826 (13) 

without abs. 

BPL 


45.647 (26) 

9.963 (13) 

11.269 (12) 

7.290 (12) 

28.960 (12) 

8.826 (11) 

19.366 (12) 

model 

PLEC 


42.999 (26) 

13.115 (13) 

21.817 (12) 

7.422 (12) 

17.094 (12) 

5.158 (11) 

17.037 (12) 

X^indf) 

LPr 


17.591 (26) 

7.245 (13) 

13.017 (12) 

4.582 (12) 

6.526 (12) 

3.181 (11) 

20.495 (12) 

with abs. 

BPLr 


25.071 (25) 

7.242 (12) 

9.870 (11) 

5.694 (11) 

4.282 (11) 

2.069 (10) 

17.474 (11) 

model 

PLECt 


43.017 (26) 

13.117 (13) 

21.823 (12) 

7.424 (12) 

17.095 (12) 

5.159 (11) 

17.042 (12) 

IN 

<1 

LP/PLt 


1.173eH-01 

3.360e-02 

7.208e-03 

2.747e-02 

4.989e-|-00 

3.267e+00 

2.331e+00 

(ndf) 

BPL/BPLt 


2.058e+01 

2.720e-|-00 

1.399e-|-00 

1.596e+00 

2.468e-|-01 

6.757e+00 

l.S92e+00 


PLEC/PLECt 


1.812e-02 

1.88Se-03 

6.614e-03 

1.464e-03 

1.646e-03 

6.277e-04 

4.025e-03 


LP/LPt 


6.166e-04 

8.546e-01 

9.323e-01 

8.684e-01 

2.551e-02 

7.070e-02 

- 

p-value 

BPL/BPLt 



9.90Se-02 

2.369e-01 

2.064e-01 


9.337e-03 

- 


PLEC/PLECt 








- 
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Table IV Fitting parameters and derived significances for the 3 flares we studied in our FSRQ sample. We mention (*) 
when a reached the lower edge of the fitting interval. 


Parameter 

Model 


3C 454.3 
(Nov-Dec 2010) 

PKS 1510-089 
(Feb-Apr 2012) 

PKS 1424-24 
(Sep 2012-Jul 2013) 


llPjunbinned 

a 

2.152 ± 0.008 

i 2.268 0.027 

2.022 ± 0.010 


LP junbinned 

p 

0.0895 ± 0.0036 

0.0451 0.0108 

0.0766 ± 0.0052 


PL 

Eb 

286 ± 8.9 

260 6.0 

940 ± 35.4 

function 

LP 

a 

2.153 ± 0.016 

2.299 0.040 

2.069 ± 0.014 

(without BLR 

LP 

p 

0.0879 ± 0.0088 

0.0569 0.0117 

0.1134 ± 0.0120 

absorption 

BPL 

Fi 

2.140 ± 0.022 

2.367 0.038 

2.010 ± 0.003 

model) 

BPL 

Fa 

2.617 ± 0.036 

2.867 0.076 

2.530 ± 0.055 


BPL 

Eb 

921 ± 0.4 

2603 912.7 

3062 ± 1.6 


PLEC 

PpLEC 

2.183 ± 0.034 

2.269 0.049 

1.903 ± 0.021 


PLEC 

Ec 

9981 ± 2111.4 

7921 2586.5 

8361 ± 1117.7 


LPr 

a 

2.157 ± 0.018 

2.299 0.022 

2.069 ± 0.014 


LPr 

p 

0.0764 ± 0.0116 

0.0569 0.0065 

0.1134 ± 0.0120 

function 

BPLr 

Pi 

2.156 ± 0.029 

- 

1.900 ± 0.019 

(with BLR 

BPLr 

r2 

2.518 ± 0.050 

- 

2.426 ± 0.036 

absorption 

BPLr 

Eb 

921 ± 0.2 

- 

1375 ± 0.2 

model) 

PLECr 

PpLEC 

2.183 ± 0.034 

- 

1.903 ± 0.021 


PLECr 

E, 

9983 ± 2113.2 

- 

8361 ± 1135.7 


LPr 


0.00703 ± 0.00551 

0.04667 0.04550 

0.00001 ± 0.00132 

a 

BPLr 


0.01059 ± 0.00591 

- 

0.00001 ± 0.00320 


PLECr 


0.00001 ± 0.00629 

- 

0.00001 ± 0.00108 


LP junbinned 


4.553 (14) 

67.828 (17) 

79.634 (15) 

X^{ndf) 

LP 


4.041 (11) 

20.807 (14) 

37.481 (12) 

without abs. 

BPL 


13.158 (10) 

15.137 (13) 

101.358 (11) 

model 

PLEC 


6.912 (10) 

37.461 (13) 

77.643 (11) 

X^indf) 

LPr 


0.928 (10) 

20.807 (13) 

37.489 (11) 

with abs. 

BPLr 


3.780 ( 9) 

- 

24.648 (10) 

model 

PLECr 


6.914 (10) 

- 

77.652 (11) 

Ax" 

LP/PLr 


3.113e+00 

6.395e-ll 

7.583e-03 

(ndf) 

BPL/BPLr 


9.378e+00 

- 

7.671e-b01 


PLEC/PLECr 


1.586e-03 

- 

9.270e-03 


LP/LPr 


7.769e-02 

- 

- 

p-value 

BPL/BPLr 


2.195e-03 

- 

- 


PLEC/PLECr 


O.OOOe+00 

- 

- 


eConf C141020.1 





























































































